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In ferromagnets, magnons may condense into a single quantum state. Analogous to supercon¬ 
ductors, this quantum state may support transport without dissipation. Recent works suggest that 
longitudinal spin transport through a thin-film ferromagnet is an example of spin superfluidity. Al¬ 
though intriguing, this tantalizing picture ignores long-range dipole interactions; we demonstrate 
that such interactions dramatically affect spin transport. In single-film ferromagnets, ’’spin super¬ 
fluidity” only exists at length scales (a few hundred nanometers in yttrium iron garnet) somewhat 
larger than the exchange length. Over longer distances, dipolar interactions destroy spin superflu¬ 
idity. Nevertheless, we predict re-emergence of spin superfluidity in tri-layer ferromagnet-normal 
metal-ferromagnet films of ~ I /xm in size. Such systems also exhibit other types of long-range spin 
transport in samples several micrometers in size. 


When matter enters a superfluid phase, it behaves like 
a fluid with zero viscosity and can support currents with¬ 
out dissipation. It has been suggested that certain ferro¬ 
magnets may exhibit spin superfluidity (SSF) [T]-[3]- The 
superfluid spin-drag properties induced by spin transfer 
and spin pumping (SP) in a normal metal-ferromagnet- 
normal metal system have recently been computed nr 
[6]. Related studies have also explored Josephson spin 
currents between magnons condensates [7]. Experimen¬ 
tal studies have suggested that the temporal decrease of 
magnon condensates is associated with SSF [8] . 

In the absence of magnetic fields, SSF is indeed an 
intriguing possibility because its realization would allow 
spin currents to propagate without significant losses over 
long distances. These spin transport properties may be 
useful for low-dissipation interconnects, spin logic de¬ 
vices, and non-volatile magnetic memory devices. Our 
work demonstrates that SSF can exist in thin-film fer¬ 
romagnetic systems, but two ferromagnets (rather than 
one) are required to cancel long-range dipole interactions. 
We do not observe signatures of long-range SSF in single- 
film ferromagnets. 

Recent works have hypothesized that easy-plane fer¬ 
romagnetic thin films exhibit SSF. In such systems, a 
monotonously precessing magnetization leads to meta¬ 
stable spin-current-carrying states whose topological 
properties protect against dissipation [3] . Spin relaxation 
induces a finite resistance proportional to the system size 
[4]. Nevertheless, ferromagnetic insulators (FIs) have ex¬ 
ceptionally low spin dissipation rates, and the spin super¬ 
current decays over a large length scale. Furthermore, the 
spin-relaxation-induced algebraic decay of the spin super¬ 
current significantly differs from the exponential decay 
of the spin current carried by spin waves [9]. Although 
magnetic anisotropy destroys the linear SSF response, 
the spin current is predicted to flow with negligible dis¬ 
sipation when the bias is sufficiently large gig. 

It is well known that long-range dipole interactions 
dramatically affect the spin-wave dispersion in thin films 
nnun]. Low-energy magnons strongly interact, and the 


coupling between them decreases algebraically as they 
spatially separate. Magnon interactions also influence 
Bose-Einstein condensation such that the condensate oc¬ 
curs at a finite wavevector around the magnon energy 
minimum [12]. Naturally, the long-range nature of dipole 
interactions can also strongly affect the SSF. However, 
for SSF to be useful, it must exist over long, hopefully 
macroscopic, length scales. 

In the previous theoretical investigations of SSF pre¬ 
sented in Refs. [m the approximate dipole field was in¬ 
cluded as an easy-plane anisotropy. However, the dipole 
interaction also has a dynamical component not included 
in Refs. IIHi It is the long-range nature of this component 
that qualitatively changes the dispersion of magnons m- 
When the system is smaller than the exchange length, the 
energy associated with the exchange stiffness dominates, 
and the system may exhibit SSF. However, the exchange 
length in FIs, such as yttrium iron garnet (YIG), is ^ 20 
nm, and dipole interactions become increasingly impor¬ 
tant at larger length scales. 

In this Letter, we investigate the complete effect of 
dipole interactions on spin transport through an FI thin 
film. We consider both square and circular devices, as 
shown in Fig. As expected, dipole interactions com¬ 
pletely alter the spin transport properties. We find that 
”SSF” can only be achieved when the system size is on 
the order of the exchange length, which implies that SSF 
is not a useful method for transporting spin information 
across sizeable distances. For example, in YIG, which is a 
widely used FI because of its low dissipation, SSF occurs 
only across distances of a few hundred nanometers. In 
comparison, typical spin-wave propagation lengths may 
reach ^ 5 mm in YIG logic devices [13]. Moreover, be¬ 
cause of dipole-induced anisotropies, a sufficiently high 
spin accumulation bias is required to induce a spin cur¬ 
rent. As with the spin resistance across the sample, this 
spin-accumulation threshold strongly depends on the ge¬ 
ometry of the system and increases with the system size. 

Nonetheless, the concept of SSF in ferromagnetic sys¬ 
tems remains useful, but not in single films, as previ- 
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FIG. 1: (color online), (a) A square FI thin film in contact 
with a spin injector (top left) and a spin detector (top right). 
A bias voltage Vi injects a spin current via a lateral spin valve 
in which the ferromagnetic leads have a perpendicular mag¬ 
netic anisotropy. SP from the FI thin film into the detector 
induces a voltage Vd- (b) The width of the contacts, ic, is 
20% of the total length of the FI sample, L. (c) Circular disk 
with diameter D. The spin injection and detection contacts 
each cover 20% of the disk area. 


ously envisaged. Instead of the single-film configuration, 
we propose a tri-layer structure. Exchange coupling be¬ 
tween two FIs via a normal metal can secure an anti¬ 
parallel configuration of the magnetizations in the two 
ferromagnets. We demonstrate that such synthetic an- 
tiferromagnets maintain long-range SSF over distances 
much greater than the exchange length. We also show 
that even when the two films differ, a spin super current 
and ultimately a long-range non-SSF spin current can 
flow over sufficiently long distances in typical realizations. 

The setup in Ref. 0] nicely illustrates SSF behavior. 
The spin Hall effect leads to spin injection. In turn, spin- 
transfer torque (STT) causes the magnetization to pre- 
cess, thereby leading to SP out of the opposite contact. 
This SP is detected via the inverse spin Hall effect. This 
geometry therefore requires the contacts to be attached 
to the thin sides of the FI. The resulting resistance per 
area can be expressed in the form of an Ohm’s law using 
the interface resistances and an internal resistance, 

^ = i/fi'f + i/flA(1) 

where and are the transverse (’’mixing”) inter¬ 
face spin conductances, and the internal spin-relaxation- 
induced resistance is = 9a/{9R9L)^ where g^ = 
2e^MsLaQ/h?^. Here, Mg is the saturation magneti¬ 
zation, 7 is the gyromagnetic ratio, is the intrinsic 
Gilbert damping coefficient, and L is the system length. 
The system exhibits SSF because the internal resistance 
Tq, vanishes when the spin is conserved {ao 0). With 
spin dissipation, the internal resistance increases alge¬ 
braically with the length of the system. 


To further utilize SSF, we suggest using a larger injec¬ 
tion area with a spin valve attached to the top of the FI; 
see Fig.[^a). Ignoring dipole interactions, Ohm’s law 0 
remains valid, but the intrinsic conductance becomes 


Qa = 


ATrMgVao 
h^Ac h ’ 


( 2 ) 


where V is the volume of the FI and Ac is the injec¬ 
tion/detection contact area. From Eqs. 0 and 0 , one 
can conclude that the SSE can be made arbitrarily long- 
range by increasing the injection area in proportion to 
the length between the detector and injector contacts. 
Without dipole interactions, the SSE is limited only by 
the contact conductances, and a spin current can flow 
over macroscopic lengths. However, as discussed below, 
dipole interactions dramatically reduce the applicability 
of this finding. 

In the geometry employed herein, a spin current is in¬ 
jected using the left contact (L), which consists of a spin 
valve with two ferromagnets, El 1 / 2 , that exhibit perpen¬ 
dicular magnetic anisotropy and are coupled to a normal 
metal, Np; see Eig. [^a). We can calculate the injected 
spin accumulation in Np, /j^l = using circuit theory 
m- Assuming an effective conductance across the NpjEI 
interface, and low spin memory loss in Np, we find 
that /ip = eVi{g^ + + ^fi)- Here, g^^) is the 

conductance of the majority (minority) electrons across 
the two Ep i/ 2 |Nl interfaces. 

This spin accumulation then drives the El dynamics 
of the local magnetization direction, m(r, t), at position 
r and time t. The spin angular momentum transported 
through the El thin film is subsequently detected by the 
right contact (R), which consists of a normal metal, Nr, 
connected to a ferromagnet, Er. The spin accumulation 
pumped into Nr is given by /J^R{r) = —hm x m|reR and 
can be measured according to the voltage, I 4 , across the 
Nr|Er junction. 

At low temperatures, SSE can be described semi- 
classically |3H6]. The magnetization dynamics are then 
described by the Landau-Lifshitz-Gilbert (LEG) equa¬ 
tion, 


ih = —7m X Heff + am x ih — a'm x m x /x/fi, (3) 


where, in the left (right) contact region, the spin accu¬ 
mulation fi = and a' = a^i^Ry^ both quantities are 

zero otherwise. The dimensionless parameter (ap(R) = 
where d is the El thickness. The lo¬ 
cal Gilbert damping coefficient is a = ao + (a', where 
a' is the spin pumping enhancement. The effective field, 
Heff, consists of the exchange field. Hex = 2 A/(M^) V^m, 
where A is the exchange constant, and the dipole field, 
H^%, which fulfills Maxwell’s equations in the magneto¬ 
static approximation, 

V X =0, V ■ + 47rM,m) = 0. (4) 
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The dipole field is related to the local magnetization 
by Green’s functions: d3r'G(r - r') m(r',t), 

where G is a 2nd-rank tensor whose elements are Gai 3 = 
-(l/47r)a2^,(l/|r-r'|) [H]. We consider an FI thin¬ 
ner than the magnetic exchange length, d < /ex = 

A/27vM‘^ ^ such that any variation of m across the 
thickness is negligible. Then, one can divide the total 
dipole field into an easy plane term. Hep = —AirMsfrizZ^ 
and the remainder of the dipole field, Hdip. 

A consequence of the dipole field H^ip is that the spin- 
wave eigenspectrum strongly depends on the spin-wave 
propagation direction relative to the magnetization m- 
At long wavelengths and no applied magnetic field, spin 
waves propagating with wave vector k\\ parallel to the 
magnetization are exchange dominated, and their fre¬ 
quency is = y/cy VSttA. In the perpendicular con¬ 

figuration, the spin waves are governed by the dipole 
interaction and uj{k±) = 4:7rMsj\/2k±d^ where k± is 
the perpendicular wave vector. SSF is associated with 
steady-state solutions of Eq. (3) where the magnetization 
has a small out-of-plane component and performs 27r pre¬ 
cessions [3]. Hence, the relative orientation of the trans¬ 
port direction and the magnetization alternates between 
the exchange-dominated and dipole-dominated regimes. 
Non-local dipole interactions are therefore of a crucial im¬ 
portance. However, the full inclusion of these interactions 
transforms the LEG equation ^ into a complicated 2^*^- 
order non-linear integro-differential equation in time and 
the in-plane coordinates, and finding its solution requires 
considerable numerical efforts. For this purpose, we per¬ 
formed graphics processing unit (GPU)-accelerated mi- 
cromagnetic simulations on several computers over a long 
time period M- 

We first consider a square YIG thin film. The square 
geometry results in two dipole-induced easy axes that 
extend diagonally across the sample in addition to the 
easy plane anisotropy. The injection and detection con¬ 
tacts cover 20% of the thin film’s surface area, as shown 
in Fig. [^b). Hence, the ratio Ac/V = l/(5d) that con¬ 
trols the internal conductance © is independent of the 
length L. We further neglect any spin-memory loss in¬ 
side the contacts because of the long spin-diffusion length 
of the Cu contacts, /^^ = 100 — 1500 nm [17]. Re¬ 
garding the Gu| YIG interface, spin-pumping experiments 
have measured transverse (mixing) conductance values 
in the range of 9cu\yig^/^^ ^ 10^^“^^ cm“^ [181 1 ^ - 
We choose = 5 • 10^^ cm“^, which, combined 

with dyiG = 5 nm, yields ~ 0.01. We also use 

An Ms = 1750 G m,A = 3.7 • 10 ^ erg/cm [20] . and 
ao = 1 ’ 10“^ [181119]. The Nr|Fr interface is assumed 
to be a tunnel interface; therefore, the SP current across 
the FI I Nr is compensated by the STT generated by fiR. 
Generally, fiR contains both ac and dc components. We 
denote the z component of the dc spin accumulation in 
the right contact averaged over the contact area by (//r) 



FIG. 2: (color online), (a) Time evolution of /Xr when 
L = 250 nm and rl — 3.0 /xeV. (b) The average dc spin 
accumulation in the right contact, (/Xr), as a function of the 
average spin accumulation in the left contact, /xl, iii samples 
with fixed lengths, L, and square geometries, (c) (/Xr) as a 
function of L. 


and investigate its behavior as a function of /xl and the 
system size. 

In the micromagnetic simulations, we start in a uni¬ 
form state and let this state evolve into a steady state; see 
Fig. I^a). Figure [^b) shows the dc spin accumulation in 
the steady state, as a function of /x^. For small val¬ 

ues oi (/Xr) is finite only above a threshold value, /x^, 
which increases with increasing system size, L. Above the 
threshold, (/Xr) increases linearly with /x^ until it satu¬ 
rates at a maximum value, /Xr^^. The vanishing (/Xr) for 
9l < dipole field that, for sample sizes 

> /ex, creates a shape anisotropy sufficiently large to pin 
the magnetic moments of the sample. For /iL ^ {ku) 
is close to the ideal value given by Eqs. 0 and Q, thus 
indicating the emergence of SSE. The saturation of (/Xr) 
at large values of kL is caused by the interplay of the 
shape anisotropy and non-local magnon-magnon interac¬ 
tions. This saturation is analogous to the leveling off of 
the cone angle in a ferromagnetic resonance experiment 
as a function of applied power due to magnon-magnon- 
interaction-mediated Suhl instabilities [2T1[22]. Note that 
because /x^^ decreases with increasing L, the interval in 
which SSE is possible, /x^ < kL < kT^^ shrinks with 
increasing L until kT^ < k^L‘ This ’’squeezing” effect re¬ 
stricts SSE to samples less than a few hundred nanome¬ 
ters in size. 

In Eig. I^c), (/Xr) is plotted as a function of system 
length, L. Eor small system sizes, (/Xr) is independent 
of L, as in the SSE Eqs. 0 and 0 . Eor intermedi¬ 
ate film sizes, the dynamic dipole field dominates the 
exchange interaction. Although spin transport is still 
possible, it is not mediated by the SSE. Eor large sam¬ 
ple sizes, dipole pinning suppresses any spin transport. 
This result demonstrates that the SSE is restricted to 
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FIG. 3: (color online). The average dc spin accumulation 
in the right contact, as a function of (a) fiL and (b) 

D for a disk-shaped thin him. as function of (c) 

and (d) D for disk-shaped tri-layer samples. The circles in 
panels (c) and (d) represent the results for symmetric tri¬ 
layer structures, where 1/2 = Mg, whereas the triangles in 
panel (d) indicate results for asymmetric layers with Ms^i /2 = 
(1 ± 0.1)Ms The solid red lines in (a) and (c) represent the 
theoretical SSF values of (/j^r). 


samples smaller than 300 nm in size and that coherent 
longitudinal spin transport over macroscopic distances in 
a single-layer sample with a square geometry is impossi¬ 
ble. 

Next, we turn to FI thin hlms with circular geome¬ 
tries. Such high-symmetry structures are chosen because 
the absence of an easy axis leads to much longer spin- 
current propagation lengths. The areas of the injection 
and detection contacts are l/5th of the total area of the 
YIG disk (see Fig. [^c)). Fig. [^a) shows that although 
spin transport is indeed possible for samples with diam¬ 
eters up to 1 /im, the SSF is still restricted to sample 
sizes < 500 nm and low values of Fig. |^b) simi¬ 
larly shows that is independent of D only for small 
values of D and jiR. In contrast with the square geome¬ 
try, no dipole pinning occurs because of the lack of easy 
axes, and long-range non-SSF spin transport may occur 
in micrometer-sized samples. 

Finally, we demonstrate that long-range SSF can be re¬ 
covered in a synthetic antiferromagnet structure. When 
two FI thin films are in contact via a thin normal 
metal, RKKY interaction can lead to antiferromagnetic 
exchange coupling between the two layers [23]. In the 
absence of an external magnetic field, the ground state 
has an antiparallel configuration with zero net magnetiza¬ 
tion. Dipolar interactions are suppressed over distances 
longer than the tri-layer thickness. Consequently, only 
the easy-plane anisotropy term survives in the thin-film 
limit. 

By applying spin accumulation to the top FI in the 
same manner as for the single-layer FI, one can induce 


a rotation of the FIs’ magnetizations that maintains 
the net magnetization close to zero. Under steady-state 
conditions, the SP current that flows from FIl into 
the spacer layer is exactly compensated by the SP 
current that flows from FI2, thereby resulting in a 
vanishing SP+STT torque, and vice versa. Hence, 
the RKKY interaction dominates the interlayer in¬ 
teraction. Writing the magnetizations as mi/ 2 (r,t) = 

(±^1 cos dz^l ”^1/2,^) , 

and assuming that {mi^2)z ^ 1 | 0 i ~ ^ 2 ! ^ the 

SSF hydrodynamic equations for the first layer are 


rhi^z 




- q;i(^i, (5a) 

[47rMs,i7 + ujE,i]mi^z + (^E,im2,z 

(5b) 


(For the second layer, interchange 1 2 .) Here, 1 ( 2 ) 

is the saturation magnetization of layer 1 ( 2 ) and ai = 
<^ 0,1 + <^ 1 , where (ao,i is the intrinsic Gilbert damping in 
FIl and a'^ = under the left (right) contact area 

and zero otherwise. The second layer is not attached 
to any external contacts, i.e., = 0 , and the damping 

is dominated by the layer’s intrinsic Gilbert damping, 
<^2 = 0 ^ 0,2 • The strength of the RKKY interaction is 
given by u;e,i( 2 ) = lJ/d 2 {i)^s, 2 {i)^ where J > 0 is the 
interlayer exchange-energy areal density [24] and di( 2 ) is 
the thickness of layer 1(2). The left and right contacts 
are attached to layer 1 and provide additional STT and 
SP, as in the single-layer cases. 

For small values of the spatial variation of mi 2 is 
small. Assuming symmetric layers, i.e. = Ms ^2 = 

Afs, 0 ^ 0,1 = (ao ,2 = and uje,i = ^e ,2 = we obtain 


(mIj) = -no. = 


9l 


9l + 9r'^ 9a 


9l, 


(6) 


niz = -fl/{2u}E + 47 rMs 7 ), (01 - ^ 2 ) ~ a 2 fl/u}E, and 
Qa is the intrinsic conductance defined in Eq. with 
the replacement ao (ao,i + <^ 0 , 2 • Equation (§1 can be 
identified as Ohm’s law Q. 

Fig. i c) shows the exact numerical result for 
as a function of for a tri-layer structure composed of 
two disks, each with thickness dyiG = 5 nm and coupling 
given by uje = 7.3* 10^^ s“^, where the dipole interaction 
is fully included. The detector signal is close to the ideal 
value given by Eq. (§, even for micrometer-sized sys¬ 
tems. Saturation also occurs in tri-layer systems but at 
much higher values of jul than for single-layer films; de¬ 
spite relatively large variations in 0, SSE remains stable 
because of screening of the dipole interactions. Eig. |^d) 
shows that spin transport is possible over much greater 
distances in tri-layer structures than in single-layer ones: 
the single-layer {jur) exhibits a 75% spin-signal reduction 
over the D interval 0.1 —2.0 /im whereas the tri-layer only 
experiences a 25% reduction over the same interval and 
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a 50% reduction at D = 4 /iui. The SSF is robust against 
small variations in the FI layer properties; see Fig. [^d). 
Our results demonstrate that tri-layer structures are able 
to support SSF currents for system sizes up to ^ 1 /im 
and long-range spin transport across samples several mi¬ 
crometers in size. 

In summary, we have investigated SSF-mediated spin 
transport in FI thin films. The dipole field qualitatively 
alters the transport properties so that single-layer SSF is 
possible only in systems less than a few hundred nanome¬ 
ters in size. Suppression of the dipole field in tri-layer 
structures enables long-range spin transport mediated by 
SSF over length scales up to ^ 1 fim and by non-SSF 
magnetization dynamics over length scales up to several 
micrometers. 
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